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Abstract 

We treat the case of a flat stress- free surface (i.e., the ground in seismological applications) 
separating air from a homogeneous, isotropic, soHd substratum overlain by a homogeneous, 
isotropic, solid layer (in contact with the ground) solicited by a SH plane body wave incident in 
the substratum. The analysis is first carried out in the frequency domain and subsequently in 
the time domain. The frequency domain response is normal in that no resonances are excited (a 
resonance is here understood to be a situation in which the response is infinite in the absence of 
dissipation) . The translation of this in the time domain is that the scattered pulse is of relatively- 
short duration. The duration of the pulse is shown to be largely governed by radiation damping 
which shows up in the imaginary parts of the complex eigcnfrequencies of the configuration. 
Three methods are elaborated for the computation of the time history and give rise to the same 
numerical solutions for a large variety of configurations of interest in the geophysical setting 
under the hypothesis of non-dissipative, dispersionless media. The method appealing to the 
complex eigenfrequency representation is shown to be the simplest and most physically-explicit 
way of obtaining the time history (under the same hypothesis). Moreover, it is particularly 
suited for the case in which modes can be excited as occurs when the incident wave is not plane 
or the boundary condition is not of the stress-free variety for all transverse coordinates on the 
ground plane. 



*Laboratoire de Mecanique et d'Acoustique, UPR 7051 du CNRS, 31 chemin Joseph Aiguier, 13009 Marseille, 
France. 
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1 Introduction 



This work is inspired by the problem of predicting the effects of earthquakes in cities. It is known 
that the most dangerous effects are produced in cities built on soft underground underlain by a 
hard substratum. A simple model of the city is considered herein in which the buildings are absent 
(i.e., the ground is flat), the soft underground is constituted by a homogeneous, soft layer overlying, 
and in welded contact with, a homogeneous, hard substratum. This configuration is solicited by 
a SH plane body wave and the object is to determine the time history of response on the ground, 
preferably in a numerically-efficient, physically-understandable manner. 

2 Space-time and space-frequency formulations 

In the following, we shall be concerned with the determination of the vectorial displacement field 
u on, and underneath, the ground in response to a seismic solicitation. In general, u is a function 
of the spatial coordinates, incarnated in the vector x and time t, so that u = u(x, t). 

We first carry out our analysis in the frequency domain, and thus search for u(x,li;), with u; 
the angular frequency. 

Fourier analysis tells us that u(x, t) and u(x, lu) are related by 



wherein it should be noted that u(x, to) is a generally-complex function, whereas u(x, t) is a real 
function. The second step will therefore deal with the computation of the integral in (^. 

3 Frequency domain analysis of the reflection of a SH plane body 
wave from a stress-free planar boundary overlying a soft layer 
underlain by a hard substratum 

3.1 Features of the problem 

Since everything is invariant with 2:3, the analysis takes place in the X1 — X2 (sagittal) plane depicted 



In this figure: Tq designates the (trace of the) interface between the substratum (half-space 
domain Oq) and the soft layer (laterally-unbounded domain il-i), and Fi designates the (trace of 
the) fiat ground. The medium above the latter is air, assumed to be the vacumn for the purpose 
of the analysis. The media in the layer and substratum are the elastic (or viscoelastic) solids 
and respectively. 

The incident plane wave propagates in toward the the interface Fq and the ground Fi. 
Since the latter is stress- free (i.e., the normal and tangential components of traction are nil on the 
boundary), the total displacement field vanishes in the region ^2 above the boundary (see fig. 
One can always choose the cartesian coordinate system so that the wavevector associated with the 
incident shear wave (the subscript T will constitute a reminder that we deal with shear=transverse 
waves in the following) lies in the xi — X2 plane. This is assumed herein and signifies that the 
displacement associated with this wave is perpendicular to the xi — X2 plane and therefore lies 




(1) 



in fig. m 
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Figure 1: Cross section view of the configuration of a stress-free fiat surface overlying a soft layer 
underlain by a hard solid substratum submitted to a SH plane wave, propagating initially in the 
substratum. 



in a horizontal plane. Thus, the incident wave is a shear wave and the associated displacement 
is horizontal; i.e., a shear-horizontal (SH) wave. Moreover, the motion associated with this wave 
is, due to the choice of the cartesian reference system, independent of the coordinate X3. This 
implies that the resultant total motion induced by this incident wave is independent of X3, i.e., 
the boundary value problem is 2D, so that it is sufficient to look for the displacement field in the 
xi — X2 plane. Actually, since we already know that the total displacement vanishes in the half 
plane above the boundary we must look for the total displacement field (hereafter designated by 
u^(x,a;)) only in Qq and Qi. 

Hereafter, we designate the density and Lame parameters in Qj by and A-', fi^ (for j = 0, 1) 
respectively. 

3.2 Governing equations 

The mathematical translation of the boundary value problem in the space-frequency domain is: 



^4/3/3(x,w) + (A'+^>; 



j = 1,2,3 , 1 = 0,1 



(5 = 1 



2, (2) 




/ = 0, 1 



(3) 




(4) 
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H^ul2{^,uj) = on Ti , (5) 

— 2(x, w) + U2 x(x, w)) = on Fi . (6) 

U2{x,uj) = ul{x,uj) on To , (7) 

^3(x,w) = M3(x,w) on To , (8) 

—Ui{x,u}) = — uj(x, w) on To , (9) 

(A° + 2/x°)m^ 2(x, w) + A%?_i(x, to) = (A^ + 2/x^)u^ 2(x, w) + A^u} ^(x, w) on Tq , (10) 

/x%3^2(x,w) = ^i^ul 2{^,uj) on To , (11) 

-M°('"i,2(x,t^) + 'U2,i(^''^)) = ~/^^(^i,2(x, w) + n2,i(x, w)) on To . (12) 

u^j^{x,uj) := tij(x, cj) — (5/o^io(^i ~ outgoing waves ; 

||x|| ^ 00 , X e fij , j = 0, 1 , / = 0, 1 , (13) 

wherein u^'^(x, u) is the (unknown) diffracted field in Qi, 6jk the Kronecker delta symbol, and: 

n3(x, uj) = Al~ exp(ik^ • x) , u\{x, oj) = uK^x, uj) = ; Vx G J^o , (14) 

k/ji = (/u^-j^ , krp2^ 5 ^2^1 ~ ; ^7^2 ^ — cos , = — , (15) 



0^ being the angle of incidence with respect to the X2 axis. 

Eq. is the space-frequency domain equation(s) of motion, (|H)- (|12|) the boundary condition(s), 
(jl^p the radiation condition, and (|14() - ()15|) the description of the incident wave. 

Until further notice, we drop the w— dependence on all field quantities and consider it to be 
implicit. 

3.2.1 Field representations incorporating the radiation condition 

As in the previous section, and on account of the outgoing wave condition(s) 1)13^ . we adopt the 
following field representations: 



/oo 
[$°+(A;i)A:i exp(fk°+ • x) + *°+(A:i)fc°+ exp(ik°+ • x)]dki , (16) 
-oo 
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/oo 
[$°+(A;i)A;°+ exp(ik°+ • x) - ^l+{ki)ki exp(ik°+ • x)]dki , (17) 
-oo 

/oo 
H°+(A;i)exp(ik°+ •x)dA:i , (18) 
-oo 



/oo 
[(^^~{ki)ki exp(ik}^" • x) + ^3~(A;i)A:^2 exp(ik|^~ • x) + 
-oo 

^'^+(/ci)A;i exp(zk^+ • x) + ^^+(A;i)4t exp(zk^+ • ^)]dki , (19) 

/oo 
[$^~(A:i)/c^2 exp(ik^" • x) - '^l~{ki)ki exp(ik;^" • x) + 
-oo 

$^+(A;i)A;[+exp(ik^+ -x) - ^+(A;i)A;i exp(zk^+ • x)] dfci , (20) 



wherein 



/oo 
[H;^"(A:i) exp(ik;^- • x) + r\~^{ki) exp(ik^+ • x)]dfci , (21) 
-oo 



kf = (^i>^it) . 4t = ±\/(^i) -(^i) , -(fci)'>0, -(A;i)^>0, 

(22) 



4^ = (fci,42) ' ^2 = ±y(4)'-(^i)' ' 3^y(4)'-(^i)'>0' »V(4)'-(^i)'>o, 

(23) 

with 

H = ^ = , , A:^ = ^ = —— . 24) 

^ 4 /aitvI ^ 4 /Z 

Note that these field representations involve nine unknown functions ^'^'^ , 'f'^^, H""*", 

H^"*". The latter will be obtained by applying the nine boundary conditions embodied 

in l|i|)-(fT2|). 

3.3 Application of the boundary condition(s) 

The use of ©, ©, ©, ®, and jn)), in EHll-dlll) gives rise to: 

$0+ = ^^0+ = $1- = ^1- = $1+ = ^1+ = ; V/ci G M . (25) 
The next step consists in using © in (|21|) to obtain: 

/oo 
[A;^~H;^~(A:i) + A;^+H;^+(A:i)] exp(iA;iXi)(iA;i = ;Vxi G M . (26) 
-oo 
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By Fourier inversion we get 

E]- =eI+ :=A^ -yheR, (27) 

whence 

/oo 
A^{ki) cos{klx2) ex.p{ikixi)dki = ;VxGr2i, (28) 
-oo 

wherein 

By the same token, we can write (|21|) as 



/oo 
A^'^{ki)exp{ik^-^X2)exp{ikiXi)dki = ;VxgOo, (30) 
-oo 

so that using (gHI) and together with (jSj) and (dU) leads to: 

;>oo 

exp(— z/c2/i) exp(iA;*xi) + / A^~^ (ki) ex.p{ik2h) exp{ikixi)dki — 

2A^{ki)cos{klh)ex.p{ikiXi)dki =0 ;Vxi G M , (31) 



oo 

oo 



/oo 
i;U°/c2A°+(A;i) exp(iA;2/j) exp{ikixi)dki + 
-oo 

/x^A;^2^^(/ci)sin(A;^/i)exp(iA;ixi)(iA;i = ;Vxi G M . (32) 



By Fourier inversion we find: 

A'- eM-ikih)S{ki - k\) + A^+ (ki) exp{ik^h) - 2A\ki) cos{klh) = ;VA:i G M , (33) 

-ifj,^kiA'- eM-ikih)d{ki - k\) + ifi^k^A^+ (ki) exp{ik^h) + fi^kl2A\ki) sm{klh^ = ;Vxi G M . 

(34) 

This system of two equations can be written as the matrix equation 

exp{ik^h) - cosik^h) \ _ ( -A'- exp{-ik\h)5{ki - k\) \ 

ifi^k^exp{ik^h) fi^klsm{klh)) \2A^) ~ \ifi^kiA'-exp{-ikih)5{ki-k\)) ' ^ ' 

It follows that: 

A^+{ki) = n°{ki)A'~5{ki - ki) , (36) 

wherein 

7z°4cOs(fcl/l) + i/ilfc2^sin(fc|/l) 1 .^^i ^ ^0);^] ^ (37) 

lj,^k2Cos{k2h) — 2/2^/^2 sin(A;2/i) J ' 

^'(^i) = ln\ki)A'"6{ki - k\) , (38) 



and 



2 

wherein 

H^k^ + fi^k 



n\ki: 



Jjfik\ cos{k\h) — i^^k\ sm{k\h) 



exp(-i4/i) . (39) 
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3.4 The scattered field 

The consequence of all this is that: 

M?'^(x, Lo) = u^'^(x, w) = ; Vx G J^o , 

/oo 
-oo 

u\{'K, uj) = uY" {"x., ^) = {yi, uj) = u\{'x., oj) = ,; Vx G , 

/oo 
^igifcixi cos(A;°X2)(iA;i = 7^l(A;i)^*-e^^l^l cos(A:^*X2) ; Vx G 
-oo 

wherein 



(40) 

(41) 
(42) 



(43) 



kl' ■= J{kif - {k\f = J{yf - (fcO sin 0^)2 , > , ^ki' > , j = 0, 1 {kl' = k\) , (44) 



7^°(^f 



li^kf cos{kl^h) + i^i^kl sm{kl'h) 



and 



lZ\k\' 



^•^^2* cos(/c2*/i) — ifi^kl^ sin(/c2*/i) 



exp[-2iA;^''/i] , 



exp[-zA:2'/i] . 



(45) 



(46) 



ji^kf cos{kl^h) - ifi^k^' sm{kl'h) 

These results, together with p4|) . indicate that the diffracted fields in Qq and f^i have the same 
(SH) polarization (vr in fig. ^ as the incident field. 

3.5 Total fields in the two media 

The preceding results show that the fields in the two media can be decomposed as follows: 



wherein 



and 

wherein 



u1{x,uj)=u^ {x,uj) + u^^{x,uj) ;VxGf^o 
^"-(x,^) = 4(x,w) = ^i-eiK^i-fc2°'^2] ^ 
ug+(x,a;) =7^°(A;i)^^-e*['=l^l+'=°'^2] ^ 

u\{x,u!) = ul~ {x,uj) + {x,uj) ;VxGr2i 



,1+ 



T 



'x •= (^i) ^4»-c»fe^i+fco"^2] 



(47) 

(48) 
(49) 
(50) 

(51) 

(52) 



Remark 

These results indicate that the diffracted field in reduces to a specularly-reflected wave u^~^ and 
the diffracted (as well as total) field in reduces to a sum of a refracted-reflected wave and 
a refracted-transmitted wave Ug^. 

Remark 

and itg"*" are plane homogeneous (body) waves so that the field in the half-space underneath 
the layer is composed of two body waves. 

Remark 

Since we consider only the case of geological interest in which the substratum is harder than the 
layer, the S-wave phase velocity in the substratum is larger than the S-wave phase velocity in the 
layer, which means that k^ < k^ and consequently /c2* is real, so that both and are body 
waves. Thus, the field in the layer is also composed of two body waves. 

Remark 

An important corollary of the previous remarks is (in the case of geological interest) that an incident 
(plane) body wave in the substratum can only excite (plane) body waves in both the substratum 
and the layer. Thus, if we want a surface wave to be excited somewhere underneath the ground, 
then we have to introduce some sort of modification of either the excitation, the nature of the 
media, or the nature of the boundary condition on the ground. 



3.6 Numerical results for the frequency domain response in the layer 

Recall that the frequency domain displacement response in the layer is of the form 

ul{x,Uj) = li3~(x, w) -|- M3^(x, w) . (53) 

with 

yi±(x,^) = 3^A^-(^^y[kl^i±kl^x2] ^ (54) 

so that one half of the total frequency domain displacement response on the ground at point 
x = (0,0) is given by 

2 " (0> 0' ^) = (^) ' (55) 

wherein: 

= dW) ' ^ ^ 

D{k\) = mcos(7u;) + 6sin(7u;) , (57) 



a = 1 , 5 = b{u}, s*) 



/i-^(cj)K" (w, 



7(0;, s*) = K^*(w,s*)-^ , 



(58) 
(59) 
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Figure 2: Spectrum of displacement response at x = (0,0). A = 1, a = 1, (3 = 4, j = 4, a = I, 
b = 0.7, corresponding to a case of separated pulses. The left-hand curves pertain to moduli, and 
the right hand curves to phases of the spectra. 



In the five figures |21 01 El El El we plot the spectrum ( A' (w)) of a Ricker pulse excitation, the 
transfer function -f^^^jy^) and the spectrum of the displacement response "3(0^1^) ^ note that 

the spectrum of displacement response contains spikes that are all the sharper and more intense 
the larger is the contrast between a and b. As will be shown further on, these sharper spikes lead 
to larger-duration response in the time domain. 



4 Time domain analysis of the reflection of a SH plane body wave 
from a stress-free planar boundary overlying a soft layer under- 
lain by a hard substratum 

4.1 Obtention of the time domain response from the frequency domain response 

We had 

/oo 
u-'^{'K,uj) ex.p{—iu}t)du! ; j = 0, 1 , (60) 
-00 

and due to the fact that u'^^{x,t) is a real function, we must have 

[n^(x,cj)]* = ul{x,-uj) , (61) 
(wherein the symbol * designates the complex conjugate operator) from which it follows that 

POO 

n3(x,t) = / ti^(x, w) exp(— iu;t)(ia; . (62) 
Jo 
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Figure 3: Spectrum of displacement response at x = (0,0). ^ = 1, a = 2, /3 = 4, 7 = 1, a = 1, 
b = 0.5, corresponding to a case of merged pulses. The left-hand curves pertain to moduli, and the 
right hand curves to phases of the spectra. 




Figure 4: Spectrum of displacement response at x = (0,0). A = 1, a = 1, (3 = 4, ^ = 4, a = 1, 
b = 0.5, corresponding to a case of separated pulses. The left-hand curves pertain to moduli, and 
the right hand curves to phases of the spectra. 
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Figure 5: Spectrum of displacement response at x = (0,0). ^ = 1, a = 2, /3 = 4, 7 = 1, a = 1, 
b = 0.1, corresponding to a case of merged pulses. The left-hand curves pertain to moduli, and the 
right hand curves to phases of the spectra. 




Figure 6: Spectrum of displacement response atx = (0,0). A = 1, a = 1, (3 = 1, ^ = 1, a = 1, 
b = 0.7, corresponding to a case of a so-called "anomalous" pulse. The left-hand curves pertain to 
moduli, and the right hand curves to phases of the spectra. 
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We shall employ (|62() or H6U|) to obtain the temporal response from the frequency response function 



U3(x, cj). More precisely, we shall be concerned with the evaluation of 



Ug (x, w) + ti3'''(x, exp{—iujt)duj ; j = 0, 1 . 



(63) 



4.2 The frequency content and time history of a of a plane transient body wave 

For a so-called upward-propagating SH-plane wave, we have 

u\{-x,uj) = , U2(x, cj) = , ul{-x,uj) = A"' (iv) exp{iklf (uj) ■ x) . (64) 

wherein 



^^rp — (^x^i ' ^2^2 ' ' — k'-p sin ^ ' — i^'j"' cos ^'j^ ^ A^^'' 



(J 



(65) 



and is the incident angle measured clockwise from the +X2 axis. In the above relations, j4*~ (lv) 
is termed the frequency- domain amplitude or spectrum of the incident SH plane wave. 
The izme history of this incident wave is 



u*(x,t) = u°~(x,t) 



u*^ {x,uj) exp{—iujt)duj 



or 



/>oo 

u*(x,t) = u°~(x,t) = 25R / u°~(x,cj)exp(- 

Jo 



-iujt)duj . 



(66) 
(67) 



4.2.1 Spectrum and time history of a Ricker pulse 

The amplitude spectrum A'^~{uj) of a Ricker pulse (Sanchez-Sesma 1985) is given by 



1 uo' 



UJ 



A^-{io) = -A 



wherein A, a and (3 are real constants (i.e., independent of uj). It follows that 



1 



so that the temporal history associated with this pulse is 



exp ( -i(]uj - — ]= A i-uj 



A- it) 



A^ (uj) exp{—iLot)duJ 



(68) 



(69) 



(70) 



More precisely: 



A'-{t) = -A 



exp 



ia;(/? - t) 



ijj 



4a2 



dijj 



-A 



4a3y^ 



up' exp 



-^{w2^4a2iu;(t-/3)} 



. (71) 
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But 
so that 



u? + \G?w{t - /3) = [(J + i2G?{lL - (5)f + 4a^(t - 0f , 



which, after the change of variables 
becomes 



1 \ 

A . ^ ^ 1 exp[-a^(t - / w^exp 



4a2 



(72) 



, (73) 



(74) 



A'-{t) = ( -^ ^^3^^ ^ exp[-a2(f _ /3)2] _ i2a^{t 



exp 



4a2 



{-Aj^-j=^[J2{t)-iW{t- f3)Ji{t)-Aa\t- pfUt)] , (75) 



wherein: 



poo 

Jo{t) = / exp 

J — oo 


4a2 


dvj 


Jo 


X) 

exp 


- ZZ72- 

4a2 




Jl(t) = 


/>oo 

/ roexp 

'— oo 


4a2 


dw = 


, 





■Mt) 



w exp 



dw = 2 / w exp 





tZ7 

'4^ 



Employing the following identities (Hodgman 1957): 

n2 







expi-^jd7? = ^ , 



7?^ exp [ —-p^ ] dr] 



we finally obtain 
so that 



or 



52;"'/-^' 

A*-(t) = ( -A ^J^ ^ [Aa^^/^-2a^/^4a^{t - P)^] exp[-a2(t - 
A'~{t) = ^ [-1 + 2a\t - pf] exp[-a2(t _ pf] . 



Remark 

The maxima of A^^{t)/A are obtained from dA^~{t)/dt = 0, i.e., 

2a^{t - P){2 - [-1 + 2a^{t - pf]] eicp[-a^{t - pf] = . 



(76) 
(77) 
(78) 

(79) 
(80) 

(81) 
(82) 
(83) 



(84) 
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Figure 7: Modulus of the spectrum function (uj) versus lo for various values of a. A = 1, (3 = 4. 



the solutions of which are 



Remark 

Thus, The minimum of A'''~{t)/A is attained at t = Iq = P and is equal to A'''~{tQ)/A = —1. This 
means that the minimum of A^~ {t)/A is independent of both a and /?. Furthermore, /3 is the instant 
at which the pulse attains its minimum, and this minimum is also the maximum of 

Remark 

A^~{t)/A = when — 1 + 2a^(t — = 0, i.e., when t = /3 ± so that 1/a is an indicator of the 
width of the main lobe of the pulse, i.e., the larger is a, the smaller is the width of the main lobe 
of the pulse. 

Remark 

The maxima of A^~{t)/A are attained at t = t_ and t = t^, and their value (0.4463) is independent 
of both a and f3. 

Remark 

The moduli of the spectra of the Ricker pulses for various values of a are depicted in fig. [Tj Note 
that these spectra do not depend on The latter only affects the phase of A'^~{io). 




(85) 



It follows that 





(86) 
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Figure 8: Time history (t) versus t for various values of a. (5 = 4: and ^ = 1. 



Remark 

The time history of the Ricker pulses for various values of a are depicted in fig. |S1 for the fixed 
value /? = 4. This confirms the fact that the larger is a, the narrower is the Ricker pulse. 



Remark 

Let G{t) be the Gaussian function 



A 



(87) 



Then 

dG{t) 



-A{t- P)exp[-a^{t- (3f] , 



(fG{t) 



A[-l + 2a^{t - /3)2] exp[-a\t - pf] = A'-{t) , 



which shows that the Ricker pulse is identical to the second time derivative of the Gaussian pulse. 



4.2.2 Spectrum and time history of a Ricker pulse plane body wave propagating in 
free space 

The wave we are concerned with is actually what was heretofore termed the " incident plane wave" . 
The plane body wave nature of the disturbance is embodied in the frequency domain function 



U3 (x,a;) = A^' {uj) eyi-p[i{k\xi - ^2X2)] = A"" (w) exp 



Xl j 



X2 Oi 



wherein 



sm( 



T ) 



1 - {giy = cOsO't 



(89) 
(90) 
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The Ricker pulse nature of the plane wave is embodied in the spectrum function 



so that the time history of this wave is given by the Fourier transform 

/oo 
A'~ (tj) ex.p[-iuJT{x,t, s')]duj 
-oo 

wherein 



(91) 



(92) 



(93) 



We assume that is non-dispersive, i.e., and Cj, do not depend on uj, so that r(x, t, s*) is 
independent of uj. Then m|(x, t) is of the same form as A^~{t) (see (jTOJ) if we replace therein t by 
r(x,t,s*), i.e., 

4(x,t) = ^^-(r(x),t,s*) , (94) 
so that employing the previous result in (|HH|) : 

n|(x, t)=A[-l + 2a2{r(x) - /3}^] exp[-a2{r(x, t, s') - f3}^] . (95) 

Remark 

The same remarks apply to this time history as to A^~(t) except that r(x, t,s*) replaces t. Thus, 
the Ricker pulse plane body wave attains its maximum (in absolute value) at r(x, t,s*) = (3, and 
the main lobe of the pulse is all the narrower the larger is a. 



4.3 Time history of the reflected and transmitted plane body wave pulses in 
the basement and layer 

4.3.1 Preliminaries 

Recall that = -j-. We encountered previously, in connection with the frequency domain response 
in Qq and Qi, plane wave functions of the type 



(x,w) := {s\uj)A'~ (uj) exp[i{k{xi ± ki'x2)] 



n^{s\uj)A'-{uj)exp 

wherein (recall that we assumed that Cj^ does not depend on uj): 
u 

r\jr\ 



lUJ I s"-^ -0 



; j = 0, 1 , (96) 



k^ 



{yy - {k\Y 



{ki/k^Y - {s^Y = v[4/4(^)]^ - 



(97) 



and, recalling that was also assumed not to depend on uj: 



n\uj,s') 



fi^K^'is') cos ( k''{uj, s*)4- ) + ifi'{uj)K^'{uj, s') sin ( k'^uj, s') 



//Ok'^*(s*) cos I K^^{uj, s*)4- ) — ifi^{uj)K^^{uj, s*) sin I ^^^(ti;, s*) 



>-jr 



(98) 
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wherein 



2/k0»(s») 



li'^K^'is') cos K^'icj, s')-^ - s') sin ^^^(a;, s')^ 



o)K-(uj,s-) ■ ( li( i\ h 



_z COS (^1^(0., + ^^^^^^^ -n (.1^(0., 



2i 

zcos(^i^(a;,50fc°^) + ""^^^^^^ sin(^i^(u;, j:) 
Thus, we are faced with the problems of evaluating the integrals 



,0+ 



(x,t)= / 

J — ( 



A^-Hi?°(w,s>xp 



(99) 



(100) 



(101) 



duj , (102) 



nj±(x,t) 



A' (uj) — — ^exp 



or 



wherein: 



/oo 
A^-(w)i2°(w,s*)exp [-iu;r°+(x,t,/)] , 
-oo 

n3^(x,t)= / A'~(u;) ^ exp [-iujT^^{^,t,uj, s')] duo 

J -oo 2 



rO+(x,t,.') := t - - k0^(/)5 + 2kO^(/)A , 

C^rjn 

ri±(x, s^) := t - T k^X.., s*)^ + k°*(/)4 • 



dio , (103) 

(104) 
(105) 

(106) 
(107) 



In the following, whenever numerical results are given, they will apply to the field on the ground 
at point X = (0, 0). 

Recall that the time domain displacement response in the layer is of the form 



so that 



ul{x, t) = u\ (x, t) + u]^{^, t) . 

!i^Ml^)=,l±(o,o,t). 



(108) 



(109) 



It is this function, related to the temporal displacement response on the ground, which will be 
depicted in the graphs that are presented hereafter. 
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4.3.2 Evaluation of M3^(x, t) for Ricker pulse excitation by a rectangle quadrature 
scheme 

The time history of the displacement field in f^i is of the form 



ni±(x,t) 



exp [—ioJT (x, t, uj, s*)] duj 



ex.p[—iujT^ (x, t, u, s*)] dcj , (110) 



wherein 



D{ijj,s^) = iacos{"fuj) +6sin(7a;) , 

a=l , b = b{u^,s') = ^ \' o- V , 7^,5 ='^ ^,^ — 
and the Ricker pulse excitation is represented by 



yl-(a;) 



Consequently 



D{—io, s*) = ia cos(7u;) — b sin(7u;) , 



(111) 
(112) 

(113) 

(114) 
(115) 



whence 



nf (x,t) 



■4 \ r°° J r exp { - r^^(x, t, cj, g')]a;} 



+ 



mcos(7u;) — 6sin(7u;) 

exp — T-^^(x, t, UJ, s^)]u;] 



iacos{juj) + 6sin(7a;) 







exp 


4a2 



duj . (116) 



This formula shows that the integrand is an exponentially-decreasing function of uj so that the 
numerical evaluation of the integral should pose no problems. In particular, we replace the upper 
limit of the integral by uJmax , where the latter is such that exp 
with the computation of 



,2 

max 



« 1, SO that we are faced 



ni±(x,t) 



iA 



4a3 



UJ 



exp — T-^^(x, t, UJ, s*)]w} 

-ia cos(7(u;, s')a;) — b{uj, s*) sin(7(ti;, s*)u;) 
exp {i[j3 — T-^^(x, t, UJ, s*)]w} 
iacos(7(a;, s'^)uj) + b{u), s*) sin(7(a;, s'')uj) 



+ 





" a;2- 


exp 


4a2 



dw . (117) 



To this end, we therefore employ the simplest method: rectangle quadrature. We divide the 
interval [0,u;ma2,'] into L equal sub-intervals, of width 5 = uJmax/I-', and centered at points uji = 
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{21 - 1)5/2 ; Z = 1, 2, L, so that 



ni±(x,t) 



iA5 



4a3 



E 



exp {-i[/3 - r^=^(x,t,u;/, s*)]a;} 



^ Liacos(7(a;i,s*)u;i) - 6(a;i,s*)sin(7(a;i,s*)u;;) 
exp - r^=^(x, t,a;;, s')]u;} 







exp 


4a2 



(itj . (118) 



This result seems to indicate that Ug (x, t) is complex. However, combining into one the two terms 
in [ ] , we find 



ni±(x,t) 



-A6 



E 



Aa^y/i^ J ^ cos^ (7(0;/, s^)oJi) + s*) sin^ (7(i^z> •5*'^) 

1^ (a + h{uji , s')) cos {7(0;^ - T^^ (x, t , , s'')loi ] + 

(a - h{u)u s')) cos {7(0;;, s*) - /3 + r^=^(x, t, wz, s*)a;/} 
which shows that u]^{ii.,t) is indeed real (at least for real 7(0;;, s*)). 







exp 





, (119) 



Eq. (|119j) is the formula we employ for the numerical evaluation of (x, t). Note that this 
formula is exact in the limits to max — > 00, L — > 00. 

4.3.3 Evaluation of n3^(x, t) for Ricker pulse excitation by a power series quadrature 
scheme 

Once again, the time history of the displacement field in Vli is of the form 



A'~{io) — ^ ' ^ exp [-iu;r^^(x,t,u;,s*)] dio 
-00 ^ 



exp [-iuJT^^{x,t,uJ,s^)] duj , (120) 



wherein 



D(uj,s^) = iacos{'juj) +6sin(7u;) , 
a=l , b = b{u,s') = ^ ' n n,v A ' 7(w, « ) = (w, s*)^ 



and the Ricker pulse excitation is represented by 



,0 ' 



A^~{u) 



:i2i) 



:i22) 



(123) 



From here on, we adopt a strategy that is different from the one in the previous section, notably 
by writing D{uj,s^) as 



(124) 
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The idea is to express D as something hke 1 — x then to employ the power series expansion 
of (1 — x)~^ to express D~^, but we have to be careful to have ||x|| < 1 in order for this series to 
converge. Thus, we write 



D{u;,s') := D-{io,s') = --{b{uj,s') - a}e'^^'^'''^^ 



v6(u;,s*) 


- 




+ a 




-a 


'b{oj,s') 





< 1 , (125) 



b{oJ, S*) + Oy 

b{uj, s*) + a 



so that: 



[D-{u;,s')]-'- 



-1 oo 



E 

m=0 
-1 oo 

E 

m=0 



b{uj, 


s'')+a\ 


b{uj, 


s^)-a) 


'b{uj 









or 



wherein 



m=0 
oo 

(6(c^,s0 + a)™ 



m=0 



ib{uj, s*) — a) 



m+l ' 



C+(^,s^) = 2 



- a 



(6(a;,s*) +a) 



m+l ■ 



< 1 , (126) 

(127) 

(128) 

(129) 
(130) 

(131) 
(132) 



At this point we recall that 



a = I , b = b[uj, s ) — • " ^ — - . " ^ 



„0 \ 2 



c^(a;) 



rO+(x,t,s*) := t - - k0^(s^)5 + 2kO^(s^)A 



■^jT.-(.,.^)J + .n.^)J 



r^±(x,t,u;,s*) :=t-s*:^TK^H^,s*)^ + '^%- ; o • 



(133) 
(134) 

(135) 
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and assume henceforth that i^i and c-^ do not depend on lo so that k = k (s*) and 



a = i , b = 6(s*) 



whence 



,m+l ■ 



Furthermore, we assume that Cq. > so that k^* is real for ah 



real for all Is* I < 1, and the condition 



gS]±^)exp[-2^7(s^u;)] 



(136) 

(137) 

(138) 
(139) 

s*| < 1. Consequently 7(0;, s*) is 
^ 1 reduces to 



It follows that the time history of the field in J7i takes the form 

°° A'- {to) 



b(s')+a 
b{s')-a 



^1- 



exp [—iuJT^^{x,t,s^)] doj 



C-{s') / ^^-(u;)exp [-M^'^(x, t, s*) + (2m + 1)7(5^)}] duj ; 

n J —00 

b{s^) + a 



m=0 



6(s*) - a 



<1, (140) 



n3^(x, t) = i 



D{u;,s^) 



exp [— iwr"*" (x, t, s*)] da; 



J]C+(s*) / A^-(u;)exp[-Mri±(x,t,s*)-(2m + l)7(s*)}]d.. ; 

n J —00 

b{s^) + a 



m=0 



b{s') - a 



or, on account of the Ricker pulse nature of the excitation, 
ni±(x,t) = 



>1, (141) 



m=0 



u;^ exp 



icj{/? - T^±(x, t, s*) - (2m + 1)7(5')} 



a; 

4a2 



dio 



bis') + a 
b{s^) - a 



< 1 , (142) 
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•uf (x,() 

(. 



A \ °^ |>00 



UJ 



2 1 



iu;{P - T^±(x, t, s') + (2m + 1)7(5*)} - ^ 



6(s0 + a 
6(s*) - a 



> 1 . (143) 



We recall here the previous result 
1 



so that 



,i± 



-A 



uP' exp 



ia;(/? - t) 



4a2 



dijj 



A [-1 + 2a2(t - (if\ exp[-a2(t - fif \ , (144) 



(x,t) = ^X;c„(.* 



m,=0 



-l + 2a2 (Ti±(x,t,s^)-/?-(s^)) 
exp f-a^ (r^±(x,t,a;,s*) -/3-(s*)) 





/6(s^) + a\ 


5 





<1, (145) 



,i± 



(x, t)=AY, C+(/) f-1 + 2a2 (r^±(x, t, s^) - /?+ (s*))' 



wherein 



m=0 



exp 



(ri±(x,t,u;,s*)-/3+(s*))^ 





/6(s*) + a\ 


1 





> 1 , (146) 

/3±(s^) := /3 ± (2m + l)7(s') . (147) 

Although these formulae are exact, they are not suitable for computation due to the presence of 
the infinite series therein. Actually, for practical (numerical) purposes, we limit the series to a 
finite (M + 1) number of terms (which is justified by the fact that the terms of the series are 
exponentially-decreasing with m) so that 



M 



m=0 



M 



(x, t)^AY^ C-{s^) -1 + 2a' (ri±(x, t, s') - (.*))' 



exp 



a' (ri±(x,t,u;,s*)-/3-(s^))' 







5 


\b{s')-aj 



< 1 , (148) 



(x,t) ^AJ2 C^i^') -1 + (ri^(x,t,s*) - /?+(s*))^ 



m=0 



exp 



a' (r^±(x,t,u;,s*)-/?+(s*))' 





/6(s*) + a\ 


5 





> 1 , (149) 



These last two formulae (which are exact in the limit M ^ 00) form the basis of what we term the 
power series quadrature method for the computation of the time history of the displacement field 
in Qi. 
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4.3.4 Evaluation of u^^{x,t) for Ricker pulse excitation by the complex frequency 
pole-residue convolution scheme 

Once again, the time history of the displacement field in is of the form 
n3=^(x,t) = / A^~{io) ^ exp [-ia;r^=^(x, t, cj, s*)] duj = 

J -oo ^ 

j^, /s exp [-itjr^^(x,t,u;,s*)] , (150) 

-oo ^[^y ^ ) 

wherein 

s*) = iacos(7Lj) + 6sin(7a;) , (151) 

„_1 h-h(:: -tM - _ ^lir, ^ 

ri±(x, t, /) := t - T ^'*(^, s*)^ + '^'''{s')^ , (152) 
and the Ricker pulse excitation is represented by 



Before going into details, we recall some general considerations. Eqs. (|15flj) - (|153|) show that the 
task is to evaluate the Fourier integral of a product of two functions F\ and F2: 

/oo 
Fi(a;)F2(a;)exp(-iu;r)da; . (154) 
-00 

We make use of the Fourier integral representations of F\{u^) and F^ioS): 

/'OO 2 /"OO 

Fi{uj) = — Fi{ti)exp{iu;ti)dh , ^2(0;) = —/ ^2(^2) exp (2(^*2)^*2 , (155) 

27r J_oo ^71" J-oo 

whose inverses are: 

/OO /"OO 

Fi{uj)exp{-iujti)dLO , F2{t2) = / ^2(0;) exp(-iu;t2)da; . (156) 
-00 J —00 

We introduce (|155l into ()154() so as to obtain: 

2 \ 2 /"OO /"OO /"OO 

- / dhFiih) dt2F2{t2) da;exp[-icj(T-ti -t2))] , (157) 

J ~oo J —00 J —00 



M(r) 



or, due to the fact that 



/OO 
da; exp [—iu>{T — ti — t2)] du = 27r(5(r — ti — t2) , (158) 
-00 
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we find 



1 r°° 

n(r) = — / Fi{ti)F2{T - h)dti , (159) 



which is a convolution integral. 

We can go a step further by expressing everything for positive times: 



1 f°° 

- [Fi{-ti)F2{T + h)+Fi{h)F2{T-h)]dh . (160) 



We now apply this formula to evaluate the time history of response. As usual, we assume that the 
media are dispersionless, so that t^^(x, t, w, s*) = T"^^(x, t, s*), whence 

ul^{^,t) = — / [Fi(-ti)F2(ri±(x,t,s^)+ti)+Fi(ti)F2(ri±(x,t,s*)-ti)](iti . (161) 

We choose 

F2{u;) = A^~{u;) , F,{co) = — ^ . (162) 

We recall the result of 

F2{t2) = A'-{t2) = ^ [-1 + 2a^{t2 - Pf] exp[-a2(t2 - Pf] . (163) 
and the Fourier inverse of Fi[ijo) is 



^i(*i)=/ TTi ^^M-i^ti)d^ ■ (164) 

J-oo i^(w,s*) 

Assuming as before that b is real, we note immediately that although the denominator of the 
integral in ()164() does not vanish for it real w, it can vanish for complex u). This suggests that the 
integral can be evaluated by use of the Cauchy theorem by appealing to a suitable integration path 
in the complex plane. Actually D vanishes at an infinite number of locations in the complex 
= Lo + uj plane, so that we prefer to proceed as following. 
In order to stress the fact that the integration variable in p64j) is real (i.e., lo ), and for other 
reasons, we re- write the integral as 

wherein 

D{uj ,s') = iacos(7a; ) +6sin(7a; ) . (166) 
We assume that D vanishes for a denumerable set of complex frequencies {tOm ; m E Z}, i.e., 

D{uj.m,s') = ; m£Z . (167) 

This suggests expanding D{uj, s') in a Taylor series around lo = uJm'- 

D{UJ, S') = D{uJm. S') + {UJ- UJm)D{uJr,,, s') + ... , (168) 
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wherein 

On account of ()167() we have 



dD{uj,s') 



du> 



whence 



1 1 



(169) 



(170) 



(171) 



Now let us turn to the issue of the actual locations of the zeros of D. We search for the complex 
roots uj of 

D(ti;, s*) = iacos(7u;) + 6sin(7a;) = iacos(7(a; + iu )) + 6sin(7(LL' + iuj ) = 0, (172) 
and assume, as was implicit (or explicitly stated), that a, h and 7 are real. Then 

D{ijO,s^) = smiy-^uj ) asinh(7u; ) + 5cosh(7C(; ) +icos(7u;) acosh(7a; )+6sinh(7u; ) 
or, owing to the fact that we have a mixture of real and complex quantitites: 



, 

(173) 



sin(7u; ) asinh(7u; ) + 6cosh(7u; 







. 



cos(7u; ) acosh(7u; ) + 6sinh(7a; ) 
Thus, we have two families of solutions, the first of which correspond to: 

sin(7Lj ) = , acosh(7u; ) + 6sinh(7a; ) = , 
and the second of which correspond to: 

cos(7u; ) = , asinh(7c<j ) + 6cosh(7c<j ) = . 
The (so-called even) solutions of the first family are: 



0J„ 



7 

27 \h + a 

and the (so-called odd) solutions of the second family are: 

,/ _ (2m + l)7r 
27 

"- = 2^^H^ 



L0„ 



m G 



(174) 
(175) 

(176) 
(177) 

(178) 
(179) 

(180) 

(181) 
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Remark 

The real parts of LOm are independent of both a and b. 
Remark 

The imaginary parts of the even frequencies are independent of m. 
Remark 

The imaginary parts of the odd frequencies are independent of m. 

Remark 

Since the ln( ) functions in these formulae are supposed to be real, the even solutions apply only 
when > and the odd solutions apply only when < 0. 



Remark 



Owing to the facts that: i) 7 > 0, ii) < |^ = 1 - ^ < 1 we have (for > 0) 



---Inf^XO 



27 v^ + " 



(182) 



Remark 



Owing to the facts that: i) 7 > 0, ii) < = 1 - ^ < 1 we have (for > 0) 



2b 



a—b 



1 , fa-b\ 
— In I r 1 < 



(183) 



27 \a + b^ 

Let us now evaluate D. We have 

-C'(wm) = 7 [-iasin(7a;m) + 6cos(7u;m)] = 7 -iasin [^{lo\^ + iu;„Jj + 6 cos [^{uj„^ + iw.^)^ 
7cos(7w^) asinh(7w^) + 6cosh(7u;^) +i7sin(7u;^) -a cosh(7w^) - 6sinh(7u;^) . (184) 

Then: 



7cos(7a;^J asinh(7u;^) + 6cosh(7u;,J^) +«7sin(7Lj^) — acosh(7u;^) — 6sinh(7a;^) 

7cos(7w^) asinh(7w^) +6cosh(7w^) = -7(-l)'" sinh(7u;^) , (185) 



7 cos(7Lj^) a sinh(7u;^ ) + h cosh(7a;°^ ) + ^7 sin(7Li;^) —a cosh(7u;^ ) — h sinh(7u;^ ) 

7sin(7u;^) -ocosh(7w^) - 6sinh(7w^) = i7(-l)"M cosh(7u;^') . (186) 
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The question that arises is whether we have to deal with either even or odd solutions. Recall that 



a = l , K^,s^)= ^o,o.(,,,.) ' 7(-,«^) = -^-,^^)^ 
and, under the (previous) assumption of a dispersionless material in the layer, 



(187) 



so that 



1\2 



-T 



i\2 



^ (/)^ [1 - is^] , 



(188) 



(189) 



or, recalling that (c^^)^ = 6^ ^ is equivalent: 



(190) 



Recall that < |s*| < 1, so that for = (i.e., normal incidence), ^ (or 6 ^ a) is equivalent 
to 

^ . (191) 
whereas for s* = ±1 (i.e., grazing incidence), b^ ^ (or 6 ^ a) is equivalent to 



(192) 



In the (geophysical) case of interest herein, i.e., dealing with a soft layer overlying a hard substra- 
tum, we have /jP > ^j} and > p^, so that we are clearly in the situation b < a for all incidence 
angles. This situation is that of odd solutions. 

Consequently, the time history of displacement in the layer is 



exp ( —iu> t] du! . 



F^{t)^iy-. : , 



Thus, we are faced with the problem of the evaluation of the integral 

1 



(193) 



/(x,t) 



■ exp {—iujt) du 



(194) 



wherein the important property to note is that iO°^ < ; Vm G Z which means that the pole of the 
integrand lies in the lower half part of the complex lo plane for all m G Z. Thus, in order to apply 
Cauchy's theorem, we consider the two contour integrals 



Ir±(x,t) = [ ex.p(—iujt)du} , 

JC± UJ-UJ^ 

wherein the contours are depicted in figs. l9land [T0l Then 



(195) 
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A Re (03) 




Iiii(oj) 



Figure 10: Contour C . 



However 



Ir+{x,t) = [ ex.p(—iujt)du} + [ exp(—iivt)dui , 

Jc+ ^-^°m JCr^- 

/^-(x, t)= / exp(— itjt)^!^; + / exp(— zci;t)(ia; . 

r 1 

I(x, t) = lim / e^T)(—iujt)dLO , 

so that, by virtue of Cauchy's theorem 

f 1 

hm /c+(x, t)= hm / exp(— ici;t)(iaj + /(x, t) = , 

ti, 

hm /^-(x, t)= hm / exp(—iu;t)(ia; + /(x, t) = —27ri Residue 



However: 



(196) 
(197) 

(198) 

(199) 
(200) 



[ ex.p{—iujt)du} < [ i? exp sin 0] 

Jet ^-^m Jo 



SO that (due to the fact that sin0 > for G [0,vr]) 

1 



hm / exp(—iujt)duj = ; t < 

R^cc Jc+ u;-UJ?n 



I{x,t) = ; t<0, 



whence (by virtue of (|199p 
Similarly 

[ ex.p{—iLvt)d(jO < f R exp [Rt sin 6] d6 

JC^ - J2n 

SO that (due to the fact that sinS < for 9 € [vr, 27r]) 



t < , R»l , (201) 

(202) 
(203) 



1 



lim 

whence (by virtue of (|200|) 

I{x,t) = — 27ri Residue 

Thus, 



exp {—iojt) duj = Q ; t > 



t > , R»l , (204) 



(205) 



= — 27ri exp(— zwt) ; t> Q . 

1° 

m 

/(x,t) = -27rzexp(-ia;^t)F(t) 



(206) 
(207) 
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wherein H is the Heaviside function {H{x) = ; X < H{x) = 1 i X > 0). 
It follows from ((TM|) that 

Due to the fact that H{-ti) = ; ti > 0, 

Fi(-ti) = ; ti > , 

whence 

1 f°° 

ul^{x,t) = — / Fi{ti)F2{T^^ix,t,s')-h)dti , 

or, more explicitly 

exp[-zu;^ii - a^(T^^(x, t, s*) - ti - 
At this point, we recall that (with obvious shorthand notation): 

< = ^i + i^^ , u;^ = (2m+l)^ , u;^' = i- hi (^^^ = a;"" < , 



(208) 

(209) 
(210) 



(211) 



(212) 



Z?(a;^,s*) = i(-l)-7^1-^cosh(7u;;^') = i(-l)™D(.*) , f)(s^) = ^^1-^ cosh (7.;°') . (213) 



a 



a 



Remark 

We notice that w^' = uj°" and D(s*) are independent of m. 
Thus, we can write (|2TT|) as 

iD Jo 

exp[-iu;^ti - a'^{T^^{x, t, s') - h - pf]dti , (214) 

wherein 



^(*i)= E (-l)'"exp(-z<ti) = exp(u;°"ti) ("ir^xp 



m=—oo 



m=—oo 



IT 



-i{2m + 1) — ti 
27 



(215) 



Thus: 
,i± 



u 



3 ,K,t) = -^ aih) [-1 + 2a2(ri±(x, t, s') -h- (if 
i\j Jo 



exp[uj°''ti-a'^{T^^{x,t,s')-ti-Pf]dti , (216) 
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wherein 

oo 



vr 



-i{2m+ 1) — ti 
27 



/ \ 00 

\ ' / m — — rv~i 



.mvr 

« (7-*i) 

, 7 



(217) 



We make use of the Poisson sum formula (Morse and Feshbach, 1953) 

exp {imdx) = — 5 f x H — 



wherein we take d = ^ and x = 7 — ti , to obtain 



(218) 



00 

J2 exp 

m=— 00 



I 7-*i 

. 7 



27 5(7-ti + 2m7) = 27 ^ 5(ii - (2m + 1)7) . (219) 

m=— 00 m=— 00 

But (recall that 7 > 0) ti = (2m + 1)7 < for m < 0, so that 

00 

27 ^ 5 (7 - ti + 2m7) ; h > , 



m=— 00 



miT 

I (7-*i) 

. 7 



(220) 



m=0 



whence 



(221) 



a(ti) =27expf-z;^ti) V5(ti- (2?n+l)7) ; > , 

V 27 

or, on account of the properties of the Dirac delta distribution, 
a{ti) = 27 ^ exp ( -i^{2m + 1)7 ) {h - (2m + 1)7) = 

m=0 V 7 / 

00 

-2^7 J]](-l)"5(ti-(2m + l)7) = ; ti > . (222) 

Consequently: 



m=0 



ni±(x,t) = f;(-l)"^ r 5{h - (2m + 1)7) [-1 + 2a2(ri±(x, t, s*) - - /3)2] x 

exp[a;°"ti -a2(T^^(x,t,s*) -ti - /?)^]dti , (223) 
or finally, on account of the sifting property of the Dirac delta distribution, 

-27^ 



ni±(x, t) = Yl (-1)" + 2a2{^'^(x, t, s^) - (2m + 1)7 - P}' 



m=0 



explu;"'' (2m + 1)7 - Q2{r^^(x, t, /) - (2m + 1)7 - . (224) 
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rectrickb 23-Apr-2006 1 1 48 alpha^l beta^l gama^l a=1 b=0.7 L=301 ommax^SO N=301 tmax=20 




I output I 

I - input I 



■ 2 4 6 8 10 12 14 16 18 20 
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Figure 11: Time history at x = (0,0) computed by the rectangular quadrature method for an 
"anomalous" pulse. ^ = 1, a = 1, /3 = 1, 7 = 1, a = 1, 6 = 0.7, LOmax = SOHz, L = 301. 



The terms of the series decrease exponentially with m so that the series can be approximated by a 
sum of M + 1 terms 

. M 

ni±(x, t) « =^ Yl (-1)™ + 2«'{ri±(x, t, - (2m + 1)7 - P}'] x 

m=0 

exp[u;°" {2m + 1)7 - a^{r^=^(x, t, s') - (2m + 1)7 - , (225) 
which is the form adopted in the numerical applications of this method. 

Remark 

As shown in the following section, the pole-residue convolution method gives rise to the correct 
solution in all cases. 

4.3.5 Comparison of the three methods for evaluating the Fourier transform inter- 
vening in the temporal response for Ricker pulse excitation 

The three methods are: i) the rectangle quadrature method, ii) the power series method, and iii) 
the complex frequency pole-residue convolution method. 

In the set of figures ^2 El El (for a so-called "anomalous" pulse), El El El (for separated 
pulses), and 1171 ITHl 1191 f for merged pulses) we exhibit the time history of displacement response at 
the location x = (0, 0) on the ground plane. 

We next choose what seems a typical geophysically-interesting situation (at least in the fre- 
quency domain) depicted in fig. 1201 The corresponding time-domain responses are given in figures 

EHEaEni 

As a last example, we choose a more idealized (and less realizable due to the very large contrast 
of physical properties it implies between the layer and the substratum) geophysical example, the 
frequency response of which is depicted in fig. 1241 The corresponding time-domain responses are 
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2 4 6 8 10 12 14 16 18 20 

l(sec) 



Figure 12: Time history at x = (0,0) computed by the power series method for an "anomalous' 
pulse. A = I, a = I, P = I, = I, a = 1, b = 0.7, M = 10. 




2 4 6 8 10 12 14 16 18 20 

l(sec) 



Figure 13: Time history at x = (0, 0) computed by the pole-residue convolution method for an 
"anomalous" pulse. ^ = 1, a = 1, /? = 1, 7 = 1, a = 1, 6 = 0.7, M = 2. 
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3ctrickb 11-Apr-2006 11 56 alpha=1 be1a=4gama=4 a=1 b=0.5 L=301 ommax=30 N=2D01 t( 

0.6 1 j-^ri 1 1 1 — - 1 
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Figure 14: Time history at x = (0, 0) computed by the rectangular quadrature method for separated 
pulses. ^ = 1, a = 1, /3 = 4, 7 = 4, a = 1, 6 = 0.5, comax = 30Hz, L = 301. 
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Figure 15: Time history at x = (0, 0) computed by the power series method for separated pulses. 
^ = 1, a = 1, /3 = 4, 7 = 4, a = 1, 6 = 0.5, M = 8. 



Figure 16: Time history at x = (0, 0) computed by the pole-residue convolution method for sepa- 
rated pulses. ^ = 1, a = 1, /3 = 4, 7 = 4, a = 1, 6 = 0.5, Af = 20. 
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Figure 17: Time history at x = (0, 0) computed by the rectangular quadrature method for merged 
pulses. A = 1, a = 2, P = A, J = I, a = 1, b = 0.1, uJmax = 'iOHz, L = 301. 




Figure 18: Time history at x = (0, 0) computed by the power series method for merged pulses. 
A = 1, a = 2, P = A, -f = I, a = 1, b = 0.1, M = 30. 




Figure 19: Time history at x = (0, 0) computed by the pole-residue convolution method for merged 
pulses. ^ = 1, a = 2, /? = 4, 7 = 1, a = 1, 6 = 0.1, M = 50. 
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Figure 20: Spectrum of displacement response at x = (0, 0). A = 1, a = 2, (3 = A, ^ = 0.5, a = 1, 
b = 0.1, corresponding to the case of a quasi-monochromatic pulse. The left-hand curves pertain 
to moduli, and the right hand curves to phases of the spectra. 



37 




10 20 30 40 50 60 



Figure 21: Time history at x = (0,0) computed by the rectangular quadrature method for a quasi 
monochromatic pulse. ^ = 1, a = 2, /3 = 4, 7 = 0.5, a = 1, 5 = 0.1, iOmax = SOHz, L = 301. 




Figure 22: Time history at x = (0, 0) computed by the power series method for a quasi monochro- 
matic pulse. ^ = 1, a = 2, /3 = 4, 7 = 0.5, a = l,b = 0.1, M = 30. 




Figure 23: Time history at x = (0, 0) computed by the pole-residue convolution method for a quasi 
monochromatic pulse. ^ = 1, a = 2, /3 = 4, 7 = 0.5, a = 1, b = 0.1, M = 50. 
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Figure 24: Spectrum of displacement response at x = (0, 0). A = I, a = 2, (3 = A, ^ = 0.4, a = 1, 
b = 0.05, corresponding to the case of a monochromatic pulse. The left-hand curves pertain to 
moduli, and the right hand curves to phases of the spectra. 
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Figure 25: Time history at x = (0, 0) computed by the rectangular quadrature method for a 
monochromatic pulse. ^ = 1, a = 2, /3 = 4, 7 = 0.4, a = 1, b = 0.05, comax = 30i/z, L = 301. 



resrickb 24-Apr-2006 13 16 alpha=2 b9ta=4 garna=0.4 a=1 b=0.05 M=150 N=1001 Imax=6 




Figure 26: Time history at x = (0, 0) computed by the power series method for a monochromatic 
pulse. ^ = 1, a = 2, /? = 4, 7 = 0.4, a=l,b = 0.05, M = 150. 
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Figure 27: Time history at x = (0, 0) computed by the pole-residue convolution method for a 
monochromatic pulse. A = 1, a = 2, (3 = 4, -y = 0.4, a = l,b = 0.05, M = 100. 
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given in figures ESI 1211 [2ZI 



Remark 

We notice that all the three methods again give the same results. 
Remark 

Note that the relatively-long duration and monochromatic nature of the temporal response in figs. 
1251 1261 1271 are due to the large Q single-spike nature of the frequency response, the latter being a 
result of the large contrast of physical properties and the fact that only one l/D spike is located 
within the significant part of the spectrum of the Ricker pulse. 

4.3.6 Discussion 

Remark 

The rectangle quadrature method, embodied in H119() . produces a purely-numerical result which 
gives no insight as to the physical nature of this response. It was proposed only as a reference 
solution by which the other two methods could be judged, at least on a numerical basis. This rect- 
angle quadrature method is certainly not optimal, even from the purely-numerical point of view, 
but obviously one of the simplest to explain and program. 

Remark 

By inspection of (|145j) and comparison with (jnH), we see that the power series method gives rise 
to an expression of the time history response to a Ricker pulse that is a sum of displaced (and 
increasingly-attenuated) Ricker pulses. This is what one would expect on an intuitive basis for a 
dispersionless configuration. Thus, it would seem that the power series method is the most ap- 
propriate one, at least in the situation in which the successive pulses are well-separated. However, 
in the case in which the successive pulses are not well-separated, intuition is lost (especially when 
a long-duration quasi-monochromatic response is produced) and the power series picture refiects 
this fact, although it still gives rise to the correct numerical response. However, the power series 
method cannot be applied when ~ as is the case in which Love modes are excited. This is the 
reason why the pole-residue convolution method was proposed. 

Remark 

The pole-residue convolution solution in ()225() expresses the time history of response as a weighted 
sum of displaced Ricker pulses (the latter would probably be distorted Ricker pulses in the pres- 
ence of dispersion). This is close to being intuitive, but what is less intuitive is the fact that the 
displacements are a function of the real part of the complex zeros of the equation D{uj,s^) = 
and the weight functions are expressed in terms of the imaginary part of the complex zeros of the 
equation D{u;, s*) = 0. 

The fact that the essential features (peak values and duration, amongst others) of the time 
history are directly-related to the complex eigenvalues of the structure is the essential result we 
were aiming at in this contribution. 

The most important parameter is the imaginary part of the eigenvalues since it regulates the 
height of the succesive Ricker pulses and therefore determines the duration of the time domain 
response. This parameter is a measure of radiation damping (which is leakage of energy into the 
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substratum, an attenuation mechanism that exists even in the absence of material dissipation in 
the layer). 

The pole-residue convolution expression of the time history appears to be similar to the one 
obtained by the power series method, but the latter method is not applicable when D{uj,s^) = 
for real eigenvalues in the absence of dissipation (i.e., the situation in which it is possible to excite 
Love modes (Groby and Wirgin 2005a,b)); moreover, the power series method does not enable one 
to predict the duration in an obvious way. 

Remark 

Some of the numerical results included in this work are rather unexpected. For instance, the time 
histories given in figs. I17 ( I18 ( ITOl have quite long durations that one would not expect to occur for 
a case in which modes cannot be excited. Actually, this long duration is due to the fact that the 
only attenuative action in this work is the one due to radiation damping. The duration would be 
shorter if material dissipation (i.e., viscoelasticity) were taken into account in the layer and/or the 
contrast between a and b were smaller. 

5 Conclusions 

The main result of this contribution is that the three methods give rise to the same solutions for a 
large variety of scattering configurations. 

The complex frequency pole-residue convolution method turns out to be the most interesting 
method since: i) it is numerically-efficient, ii) it is explicit as concerns the understanding and 
quantification of the duration of the time domain response, iii) it can be employed even in the case 
in which genuine resonances (due to mode excitation) are produced. 

The part of this study concerning the complex frequency pole residue convolution method 
constitutes a correction of its counterpart in our previous publications (Groby and Wirgin, 2005a) 
and (Groby and Wirgin 2005b). A somewhat similar approach, although applied to a fluid layer in 
a fluid host, is that of (Conoir, 1987). 

Work remains to be done to take into account dispersion and damping of the material in the 
layer. 

The natural follow-up of this study is to elucidate theoretically the nature of the time histories 
of response not only for the case (the one treated herein) in which the configuration is unable to 
excite (e.g.. Love) modes, but also in the case in which such modes can be excited (Groby and 
Wirgin 2005a,b). 
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